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Abstract 

Several classes of self-similar, spherically symmetric solutions of rela- 
tivistic wave equation with a nonlinear term of the form sign(<^) are pre- 
sented. They are constructed from cubic polynomials in the scale invariant 
variable t/r. One class of solutions describes a process of wiping out the 
initial field, another an accumulation of field energy in a finite and growing 
region of space. 
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1 Introduction 



Symmetry transformations for a given equation act within the space of its solu- 
tions. Especially interesting are fixed points of such transformations, i.e., the solu- 
tions that remain invariant. In the case of continuous symmetries they depend on 
a reduced number of independent variables. This significantly simplifies the task 
of finding the solutions, the task that so often is formidably difficult in the case 
of nonlinear field equations. Apart from the rudimentary examples of rotational 
and translational symmetries, in certain models there is also a one-parameter scal- 
ing symmetry. The corresponding invariant solutions, called the self-similar ones, 
play an important role in mathematical analysis of nonlinear evolution equations, 
as well as in physics oriented investigations, see, e.g., the books (T), and a 
sample of recent research papers 0. 

The present paper is a sequel to papers flU, where self-similar solutions in the 
signum-Gordon model with a real scalar field in 1+1 dimensions were analyzed. 
Rich variety of such solutions was found. Certain related models were investi- 
gated in BH. Derivation and a general discussion of the signum-Gordon model 
can be found in @J. The model is very interesting, there is no doubt that its in- 
vestigations should be continued in several directions, such as finding classical 
solutions of various kinds, applications, or quantization. 

As is well-known, field theories in 1+1 dimensions have rather special prop- 
erties. Therefore, it is not clear whether the amazing picture of the space of self- 
similar solutions found in flU is valid also in higher dimensions. This is our main 
motivation for the present work. Another attraction is the possibility of obtaining 
several exact analytic solutions to the signum-Gordon equation in 3+1 dimen- 
sions, and to have new insights into rather intriguing dynamics of the scalar field 
in that model. 

We have found several families of self- similar, spherically symmetric analytic 
solutions. They are composed from cubic polynomials in the scale invariant vari- 
able u = t/r. Our results show that there is a qualitative similarity between the 
3+1 and 1+1 dimensional cases, but significant differences appear. First, the lack 
of translational invariance in the radial variable r results in the absence of counter- 
parts of certain classes of 1+1 dimensional solutions. Second difference, perhaps 
not so unexpected, is that the pertinent calculations are a bit more complicated 
because we have to look for zeros of cubic polynomials, while in the 1+1 dimen- 
sional case only the second order polynomials were present. The exact solutions 
we have found provide examples of non-trivial evolution of the scalar field <p. Es- 
pecially interesting seem to be solutions of the types II and III, presented in 
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Section 3, which show how the field settles exactly at the vacuum value ip = 0. 

The plan of our paper is as follows. In Section 2 we describe the method 
of constructing the self-similar solutions in the signum-Gordon model. Explicit 
solutions in the 3+1 dimensional case are presented in Section 3. Section 4 is 
devoted to discussion of our results. 



2 Scale invariant Ansatz and the method of construct- 
ing the solutions 

The signum-Gordon equation for the real scalar field (p in D + 1 dimensional 
space-time reads 

d^cp(x) + sign(^(x)) = 0. (1) 

The sign function takes the values ±1, 0, sign(0) = 0. It is clear that Eq. (1) im- 
plies that the second derivatives of (p may not exist if ip vanishes at isolated points, 
but even at such points one-sided second order derivatives do exist. The proper 
mathematical framework for discussing equations of such kind is well-known: 
one should consider the so called weak solutions, see, e.g., H. All solutions 
constructed below have been checked in this respect. 

From a given solution (p(x) of Eq. (1) one can obtain one-parameter family of 
solutions of the form 

<p x (x) = AV(^), (2) 

where A > is a constant. The solution is self-similar if (f\(x) = ip(x) for all 
A > 0. Taking A = r, where r is the radius, and assuming the spherical symmetry, 
we may then write 

V (x) = r 2 G(-), (3) 
r 

where t = x° is the time. With this Ansatz, Eq. (1) gives 

{u 2 - l)G" — (D + l)uG' + 2DG - sign(G) = 0, (4) 
where u = t/r and G' = dG/du. This equation has the particular solutions 

which are trivial, but nevertheless play important role below. 
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Let g(u) be a solution of the auxiliary linear (!) equation 
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l)g"-(D + l)ug' + 2Dg = 0. (5) 



Then G+(u) = g(u) + 1/(2D) is a solution of (4) on the interval of u defined by 
the condition G+(u) > 0, and G-(u) = g(u) — l/(2D) on the interval in which 
G-(u) < 0. It turns out that patching together a number of such partial solutions 
and, in some cases, including also the trivial solution G — 0, one can cover the 
whole interval < u < oo. In this way we obtain a self-similar solution of Eq. 
(1) valid for all r e [0, oo) and t G [0, oo). The values (7(0), G'(0) determine the 
initial data for ip(t, r): 

p(0,r) = r 2 G(0), (%>)(0,r) = rG'(0). (6) 

When patching the partial solutions we demand continuity of and also con- 
tinuity of G'(u), unless u — 1. G"(m) at u = 1 does not have to be continuous 
because G" in Eq. (1) is multiplied by the factor u 2 — 1 which vanishes at that 
point. This is not surprising because u — 1 corresponds to the light-cone r = t, 
the characteristic hypers urf ace for equation (1). 

Equation (5) has two linearly independent solutions: 

oo 

9l (u) = Du 2 + l, g 2 (u) = Y J CiU 2l+ \ (7) 

1=0 

where the coefficients q are determined from the recurrence relation 

_ (2<-l)(2f + l-£>) 
Q+1 " 2(/ + l)(2/ + 3) (8) 

Note that in the case of odd space dimension D the series for g 2 is in fact a poly- 
nomial of order D, because q = for all I > (D — l)/2. In particular, g 2 — u 
for D — 1, (?2(w) = m 3 + 3-u when D = 3, and g 2 = u 5 — 10-u 3 — 15-u for D = 5, 
apart from overall multiplicative constants. 

In the next Section we consider in detail the D = 3 case. 



3 Explicit self-similar solutions in the case D = 3 

In the case of three space dimensions the general form of the partial solutions 
reads 

G = ±- + a(3u 2 + l) + f3(u 3 + 3u), (9) 
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or equivalently, 

G = ±l+ 1 (u-l) 3 + 5(u + l) 3 , (10) 
6 

where a, /3, 7, 5 are constants. The sign in front of 1/6 has to be equal to sign(G). 

Let us begin our exploration of the space of the self-similar solutions by check- 
ing their asymptotic behavior at u — > 00. It turns out that G has to have a constant 
sign if values of u are large enough - we prove in the Appendix the following 
lemma. 

Lemma. Suppose that u\ > 1 is a real zero of the polynomial (9). Then, that 
polynomial does not have any real zeros larger than u\. 

It follows that in the region u > u\ the solution has the fixed form (9) with certain 
fixed coefficients a, (3. The just excluded option was that of having an infinite 
sequence of polynomials with alternating signs. Now, let us recall that the region 
u — y 00 corresponds to r — > (t > 0). Therefore, the polynomial u 3 + 3u would 
give a singular field p = r 2 G ~ 1/r with a divergent, non integrable at r = en- 
ergy density. For this reason, we put j3 = 0. Note also that by taking into account 
the symmetry p — > —<p, we may consider just two cases of the asymptotic form 
of G at large values of u: G(u) > 0, or G(u) = 0. The latter case is considered 
in the second and third subsections (solutions of the types II and III). The former 
case is named type I. 

3.1 Solutions of type / 

In this case, for arbitrary large values of u 

G(u) = 1 + a(3u 2 + 1) > 0, 
6 

as pointed out in the preceding paragraph. This implies that a > 0, and then such 
G{u) does not have any zeroes down to u — 1, where it can be glued with another 
polynomial of the form (9) or (10). We conclude that all type I solutions on the 
interval u e [1, 00) have the form 

G 00 (u) = ^ + a 00 (3u 2 + l) (11) 

with an arbitrary constant > 0. 

At the point u — 1 the solution G^ can be glued with another polynomial (9) 
or (10), taken with the +1/6 because G^l) = 1/6 + > 0. Let us denote 
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such polynomial by G + and its constants by 7+, S + (here we prefer the form (10)). 
The gluing condition at u = 1, Goo(l) = G + (l) gives 5 + = a^jl. Hence, 

G+(u) = \ + 1+ (u - l) 3 + \ aoo (u + l) 3 . (12) 
6 2 

Now we have to determine the lower end of the interval on which G + is the solu- 
tion (the upper end is u = 1). The first possibility, denoted as la, is that G + > 
for all u E [0, 1). This occurs if 7+ < 1/6 + aoo/2. In this case G + and G^ cover 
the whole interval [0, 00), and these functions together give the complete solution. 
It has the shape sketched in Fig. 1 with the dashed line. 




Figure 1 : Solutions of the type la (the dashed line) and lb (the solid line) 



The other possibility, denoted as lb, is that there exists u £ (0, 1) at which 
G + (u ) = 0. Simple calculation gives in this case 

7+ = 771 r 3 • (13) 

6(1 - u y 
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Moreover, it turns out that u has to be the first order zero of G + . Therefore, at 
the point u we may glue G + with a negative polynomial 

G_ = -i + 7 „( M -l) 3 + 5_( M +l) 3 <0, (14) 
o 

and not with the trivial solution G — 0. The matching conditions 

G_(«o) = 0, G , _(uo)=G'4u ) 

give 

U + 3^00(1 + M ) 3 r «oo 1 MC s 

6(1 -u ) 6 2 6(1 + u ) 2 

It turns out that G_ (iz) found above does not have any zeros in the interval 
[0, uo). Therefore, the three functions: for u > 1, G + for u < u < 1, and 
G_ for < u < uo, form the complete solution of Eq. (4). It is sketched in Fig. 1 
with the solid line. 

The values of G(0), G'(0), which specify the initial values of the scalar field 
through formula (6), read as follows. 
Type la: 

G + (0) = l-l + + ^, G'_(0) = 3 7+ + ^. 



Type lb: 
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G_(0) = ~-7- + *-, G'_(0) = 3( 7 _ + (J-). 
o 



Varying in the interval [0, 00) and u in the interval [0, 1) we obtain certain 
sets in the (G(0), G"(0)) plane. They are shown in Fig. 2. We have denoted them 
la, lb, identically as the corresponding solutions. 

The regions —la, —lb are obtained by the reflection in the origin, (G(0), G'(0)) — >■ 
(— G(0), —G'(0)), related to the symmetry ip — > —ip. The half-infinite curve that 
separates the region lb from lib is obtained for = 0. It has the following 
parametric form with u G [0, 1) as the parameter: 

G(0) = ^A^-^^-l G{0) = *^+ U ° 



6(l + « ) 2 6(1-m ) 3 6' w 2(1 + u ) 2 2(1 -Mo) 3 ' 
The curve starts at the point (0, 1/2) that corresponds to u = 0, and it approaches 
its asymptote, that is the straight half-line 



G(0)(a) = ~-a, G'(0)(<7) = ^ + 3<7, <tG[0,oo), 



when u — > 1 . 
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Figure 2: The map of the self-similar solutions. The regions la, I la are separated 
from lb, lib, respectively, by the G'(0) axis. The infinite region 77/ extends 
over all area encompassed by the stripes I la, lib, —I la, —lib, except for the 
origin (0, 0) which corresponds to the trivial solution G = 0. The stripe lib 
asymptotically approaches the stripe —Ila, similarly Ila approaches —lib. The 
meaning of the dashed straight line is explained in Section 4 
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3.2 Solutions of type II 

Solutions of the type II are obtained by taking the trivial solution G(u) = for 
u > Ui > 1, and gluing it at the point ui with the polynomial (10) taken with the 
+ sign (the other sign can be obtained by the symmetry transformation ip — > —if). 
Because u\ ^ 1, the matching condition includes the first derivative, 

G( Ul ) = 0, G'(ui) = 0. 

Solving these conditions we obtain the following polynomial 

, {ui - u) 2 (2 Ul u + uj-3) 

G d{U) = — "a TTo • ( 16 ) 

It has strictly positive values in the whole interval u G [l,ui). At u — 1 that 
polynomial is glued with another positive polynomial G+(u) of the form (9). The 
matching condition G d (l) = G + (l) gives 

G+( .) = i-M«-i)'-g^, (17) 

where a + is an arbitrary constant. Similarly as for the type I solutions, there are 
two possibilities. The first, denoted as 1 1 a, with G + (u) strictly positive on the 
whole interval (0, 1], occurs when a + > —1/6. In this case we already have the 
complete solution: G = for u > ui, Gd for u £ [u\, 1], and G + given by 
formula (17) for u e [0, 1]. It is sketched in Fig. 3 with the dashed line. The 
values of G + , G' + at u — 0, i.e., 

G + (0) = l + a + , G" + (0) = -3a H ' 



6 ' ^ ~ + ( " J 2(1 + Ml ) 2 ' 

where a + e [— 1/6, oo), u\ G (1, oo), fill a semi-infinite straight- linear stripe in 
the (G(0), G"(0)) plane. It is denoted as Ila in Fig. 2. 

The second possibility is that G+(u) has a zero at some u G (0, 1). If this is 
the case, then 

= u (u 2 + 3) 1 

a+ 6(l- Mo ) 3 (l+«i) 2 6(1 - Mo ) 3 ' 
Formula (17) implies that the zero can only be of the first order. The matching 
conditions at u determine the two coefficients in the negative polynomial 

G-{u) = -i+ 7 _(M-l) 3 + Ml + M ) 3 (18) 
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which is glued at that point with G + (u). Simple calculations give 

up {1 + UqY 1 1 

7 " 6(1 -u ) 3 12(l + Ml ) 2 (l-n ) 3 ' " 6(1 + n ) 2 12(1 + Ul f 

It turns out that G_(u) is strictly negative in the whole interval u G [0, w ) 
for all choices of u E (0, 1), U\ E (1, oo). Thus, we have obtained again the 
complete solution: G = for u > u\, Gd for u E [1, ui], G + for u E [uo, 1], and 
G_ for m G [0, ui\. We denote it as lib. It is sketched in Fig. 3 with the solid line. 




G(u) 



-0,2 



-0,3 



-0,4 



Figure 3: Solutions of the type I la (the dashed line) and lib (the solid line) 



The points (G_(0),G'_(0)) with u E (0,1), Mi > 1 fill in the (G(0),G'(0)) 
plane the region denoted as lib, see Fig. 2. It borders the region I la along the 
segment (3/8, 1/2) of the G'(0) axis. The region lib is bounded from below 
by the curve obtained by putting u\ — 1 in formulas for G_(0), G'_(0) for the 
solution (18). The curve is obtained in the parametric form, namely 

u (2 + u ) _ Uq(1 + u ) = 1 3u -l 

" U 6(l + Mo ) 2 24(l- Mo ) 2 ' " U 2(l + % ) 2 8(l- Mo ) 2 ' 
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where u varies in the interval (0, 1). The asymptote of this curve, approached 
when u —s- 1—, coincides with the upper boundary of the region —I la in Fig. 
2 (the semi-infinite straight line that starts from the point (0, —3/8) on the G'(0) 
axis). 

3.3 Solutions of type III 

Thus far we have not found solutions for which the points ((2(0), G'(0)) would lie 
in the central part of the map presented in Fig. 2. In fact, we have not exhausted yet 
all possibilities for having the solutions of the type II because we have assumed 
that u\ > 1. Simple calculations show that the partial solution of the form Gd, 
formula (16), does not exist if u\ < 1. Thus, we are left with the last possible 
choice, namely u\ — 1. In this case G(u) = for u > 1. The matching condition 
with the trivial solution G = at u = 1 does not involve the first derivative. For 
his reason, we regard such solutions as essentially different from the ones already 
constructed, and we classify them as the type III. 

As for the general form of such solutions we are guided by the results obtained 
for the signum-Gordon model in 1+1 dimensions, [4] . Thus, in the present case we 
expect an infinite sequence {Gi, G 2 , G 3 , . . .} of polynomials of the form (10) with 
alternating signs, glued together at the points u , w 1; u 2 , ■ ■ ■ lying in the interval 
(0, 1) and converging at u = 1. These points are the first order zeros of the 
pertinent polynomials, and the matching conditions include the first derivatives. 
Unfortunately, we have not been able to fully construct these solutions in exact, 
analytic form because the cubic polynomials pose significant problems in explicit 
calculations. One can relatively easily construct several first polynomials with the 
desired properties, but the infinite sequence of cubic polynomials is a different 
matter. 

Numerical calculations support our expectations. Example of such a numeri- 
cal solution is shown in Fig. 4. It has been calculated up to the point u = 0.9996, 
determined by numerical accuracy of the computation. Of course, the numerical 
results for this kind of a problem, where we need to check an infinite sequence of 
polynomials, are not decisive, even if they are quite suggestive. For this reason, 
the given above description of the solutions of the type III has, strictly speaking, 
the status of a conjecture, as opposed to the status of the solutions of the types I 
and II. 
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Figure 4: Example of numerical solution of the type III 



4 Discussion 

1 . We have found the three families of self- similar solutions of the signum-Gordon 
equation in the D = 3 case. Comparing with the analogous results in 1+1 dimen- 
sions flU, we see certain similarity. In particular, in both cases there exist solutions 
of the type III (i.e., solutions vanishing for all u > 1) and of the type 77 (vanish- 
ing for all u > ui > 1). This seems to be a universal feature of the signum-Gordon 
equation. Solutions of the type / were not found in flU, but this may be a conse- 
quence of certain ad hoc restrictions on the form of solutions adopted there. On the 
other hand, in the D = 3 case we have not found spherically symmetric solutions 
that would correspond to the presented in |4j solutions with a negative velocity v. 
This can be related to the lack of translational invariance in the radial coordinate r. 



2. The solutions presented above describe quite interesting processes, especially 
the type II solutions. The region u > u\ > 1 in which ip = corresponds 
to r < t/u\. Thus, the field retreats completely from the region surrounding 
the origin. The radius of the spherical region with the vanishing field increases 
linearly with time. The corresponding radial velocity is equal to l/u\, and it is 
smaller than 1 . The presence of such subluminal velocities shows that actually 
the model should not be regarded as massless, in spite of the presence of the 
scale invariance. In truly massless models, such as the free massless scalar field 
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discussed in the next point, wave fronts always move with the velocity ±1. 

The type I solutions give (p(r, t) = (1/6 + a^r 2 + 3a 0O t 2 in the region r < t. 
Thus, in this case we have an accumulation of the field energy around the origin. 

3. It is interesting to compare the solutions described in Section 3 above with self- 
similar radial solutions of the free wave equation which is obtained by dropping 
the signy? term from Eq. (1). It is clear that the general solution of the free wave 
equation is given by linear combination of the functions gi(u), gi(u) introduced 
in Section 2. It turns out that solutions of the type II do not exist. Solutions of 
the type III have the form G(u) = a(l — u) 3 for u < 1, G(u) = for u > 1, 
where a is an arbitrary constant. On the plane ((2(0), G'(0)) these solutions are 
represented by the points lying on the straight line G'(0) = — 3G(0), shown as 
the dashed line in Fig. 2. All other points in that plane would correspond in the 
case of the free wave equation to solutions of the type I, which have the following 
form: G(u) = 2(3(3u 2 + 1) for u > 1 and G(u) = (3(1 + u) 3 + a{l - uf for 
u < 1, where j3 ^ 0. Thus, the rather nontrivial shape of the region occupied 
by the solutions of the type III in Fig. 2, as well as the presence of the regions 
I la, lib, reflect the presence of the sign(y?) term. 



5 Appendix. The proof of Lemma 

1. If ui is the second order zero, the polynomial has the form (16) modulo the 
overall sign. The factor 2uiU + u\ — 3 is strictly positive for all u > ui because 
ui > 1. Therefore Gd does not vanish for any u > u\. 

2. In the case ui is the first order zero, we have G(ui) = 0, G'(ui) ^ 0. Because 
of the symmetry if — > —if we may consider only the case G'(ui) > 0. Then, for 
u > ui we have the polynomial 



G+(u) = ---}- 
+y ' 6 6(1 



3u< 



3«?) 



+ (3 + (u 3 + 3«i) 



u 3 + 3u 1 + 3m 2 



u{ + 3«i 1 + 3uj 



where 



/?+> 



Ui 



3{u 2 - If 



in order to ensure that G' + (ui) > 0. Direct calculation shows that G'[(u) > for 
all u > ui. Therefore G' + (u) > for all u > ui, and in consequence the values 
of G + (u) do not return to in the region u > u\. 
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